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FOREWORD 


This  is  the  twelfth  report  in  the  Data  Reduction  Program  Documentation 
series.  It  is  date  recording  to  the  date  of  completion  of  the  documentation.  No 
implication  is  made  that  this  program  will  not  subsequently  be  modified,  amended, 
or  superseded;  on  the  contrary,  the  history  of  radar  data  processing  is  one  of  con¬ 
tinuous  evolution  of  techniques,  and  it  is  unrealistic  to  assume  that  steady -state  has 
been  reached. 

The  preparation  of  reports  in  this  series  is  under  the  Editorship  of  Charles  R. 
Berndtson  of  Lincoln,  anu  of  D.  Nessman  and  R.  French  of  PhUco-Ford  Corporation. 
Inquiries,  suggestions,  corrections,  criticisms,  and  requests  for  additional  copies 
should  be  directed  to  C.  R.  Berndtson. 

The  principal  contributor  to  this  report  was  L.  C.  'Zemke  (Philco-Ford). 

Due  to  the  intricate,  evolutionary  manner  in  which  the  programs  came  into  being, 
the  editors  regret  that  it  is  in  general  impossible  to  give  due  credit  to  all  --  mathema¬ 
ticians  or  radar  analysts  or  programmers  --  who  contributed  to  the  definition  and 
writing  of  the  programs. 


Alan  A.  Gromeystein 
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COMMON  SYMBOLS  AND  ABBREVIATIONS 


(The  units  given  for  certain  quantities  are  the  units  commonly  used  for  those 
quantities,  unless  otherwise  noted. ) 


ADT 

ALCOR 

ALTAIR 

Alt 

APS 

ARS 

ARTP 

ATC 

Avg 

Az 


CADJ 
C -band 

DBLT 

DCO 

El 

EOF 

GMT 

h 

Hz 

IF 

in 

IRV 

LC 

lsb 

min 

NB 

NRTPOD 

POD 

Phase 

PRF 

PRI 

pps 

pts 


ALCOR  Data  Tape 

ARPA -Lincoln  C-band  Observables  Radar 

ARPA  Long-Range  Tracking  and  Instrumentation  Radar 

Altitude  (km) 

Average  Pulse  Shape 
ALTAIR.  Recording  System 
ALTAIR  Real  Time  Program 
Angle  Track  Console 
Average,  Averaging 
Azimuth  (deg) 

Speed  of  Light 

Adjusted  Calibration  Constant  (db) 

ALCOR  frequency,  5664  MHz  (NB)  and  5667  MHz  (WB) 

Wide  Band  Pulse  Doublet 

Designations  and  Communications  Operator 

Elevation  (deg) 

End  of  File 

Greenwich  Mean  Time 

Hours 

Hertz 

Intermediate  Frequency 
Inches 

Inter -Range  Vector 

Left  Circular  Polarization 
Least  Significant  Bit 

Minutes 

Narrow  Band 

Non -real  Time  Precision  Orbit  Determination  Program 

Project  PRESS  Operation  and  Data  Summary  Report 
Presented  in  deg 

Pulse  Repetition  Frequency  (pps) 

Pulse  Repetition  Interval  (s) 

Pulses  per  second 
Points 

vii 


R 

• 

R 

rad 

RC 

RCS 

RF 

RGC 

RTC 


Range  (km) 

Range  Rate  (km/s) 

Radians 

Right  Circular  Polarization 
Radar  Cross  Section  (dbsm) 
Radio  Frequency 
Receiver  Gain  Control 
Range  Track  Console 


s  Seconds 

SD  Standard  Deviation  of  Wake  Velocity 

w 

SDBLT  Wide  Band  Slaved  Pulse  Doublet 

S/N  Signal -to -noise  Ratio 

T  Time 

TAL  Time  After  Launch  (s) 

TGC  Transmitter  Gain  Control 

Tr  Traverse  Angle  (deg) 


UHF  ALTAIR  Frequency;  415  MHz 

V  Velocity 

Vd  Doppler  Velocity 

V  Mean  Wake  Velocity 
w 

VHF  ALTAIR  Frequency;  155.5  MHz 


WB  Wide  Band 

WBS  Wide  Band  Slaved 

WTR  Western  Test  Range 

8  Total  Off-axis  Angle  (deg) 

A  Wavelength 

*  Denotes  Multiplication 


FLOW  DIAGRAM  SYMBOLS 


PROCESS,  ANNOTATION 


CD 


TERMINATOR 


SUBROUTINE:  where  NAME  I*  the  entry 
call  into  the  luhroutine 


CONNECTOR:  where  P  specifies  a  page  in  the 

flow  diagram,  and  L  designates 
a  statement  number  in  the  program 
listing  or  a  reference  point  in  the 
flow  diagram 

CONNECTOR:  where  X  implies  a  continuation 
of  the  diagram  to  the  neat  page 


INPUT/OUTPUT  OPERATION 


MAGNETIC  TAPE 


PUNCHED  CARD 


DISK 


ix 


ALTCEP 


DESCRIPTION 


ALTCEP  examines  an  ALTAIR  catalog  tape  for  radar  status,  performance, 
and  target  fLCS.  It  produces  Voting  (at  a  nominal  interval  of  0.1  s)  of  status  and 
RCS  data  and  plots  of  RCS  v^TAL  (VHF  and  UHF  LC  and  RC). 


The  computation  of  metric  data  depends  on  the  ARTP  version  used 
at  Kwajalein.  Since  12  March  1971,  R,  Am  and  El  are  corrected  for  known  bias  errors 
prior  to  recording  on  the  catalog  tape.  Before  12  March  1971,  R,  Az,  and  El  were 
corrected  for  bias  errors  by  the  ALTCEP  program. 

M 

Vital  format  tables"  are  tested  through  a  call  to  CHEAD,  and  if  any  are 
missing  the  program  terminates. 

The  following  paragraphs  explain  each  item  listed  by  ALTCEP. 


GMT 


ALTCEP  uses  GMT  from  two  locations  on  the  catalog  tape. 

The  GMT  found  in  FMRDRM,  Items  1  and  2,  is  used  with  the  metric 
listing.  GMT  found  in  FMCATF,  Items  3  and  4,  is  used  with  the  ARS  target  listing. 

GMT  (FMCATF)  is  compared  to  GMT  (FMRDRM),  and  if  they  differ 
by  >  1  s,  the  record  is  skipped  and  an  error  message  printed. 

B.  TAL 


TAL  is  GMT  -  lift-off  time.  Lift-off  time  may  be  input  on  data  cards 
or  read  from  the  tape.  If  lift-off  time  is  not  input  on  data  cards,  ALTCEP  examines 
Calibration  Record  FMGLOT,  Item  1,  to  determine  if  the  lift-off  time  in  FMGLOT, 
Items  2  and  3,  is  valid.  If  lift-off  time  is  valid,  it  is  used.  If  it  is  not  valid,  it  is 
set  equal  to  zero,  and  TAL  equals  total  GMT  s. 


^ Found  in  FMHDRD,  Items  10  and  11. 

## FMCATF,  FMCIDB,  FMCTDB,  FMRDRD,  and  FMRDRM. 
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C.  Alt 


2  2  i 

Alt  =  (R  +  R  +  2RR  sin  El)^  -  R  ,  where  R  is  earth  radius, 
e  e  e  e 


D.  R 


1.  Prior  to  12  March  1971 


R  =  fIR(85)+IR(86)  +  I11(3)]  -AR 


where 


!r  (85)  +  IR  (86)  is  the  unambiguous  range  found  in  FMRDRM, 
Items  85  and  86 

I11  (1 2 3)  is  range  bias  found  in  Calibration  Record  FMRR11, 
Item  3 

AR  is  tropospheric  refraction  correction 
2.  On  or  After  12  March  1971 
R  =  f  IR  (85)  +  IR  (86)1  -AR 


R  is  found  in  FMRDRM,  Item  84. 


1.  Prior  to  12  March  1971 


Az=Ir(13)+Ic  (3) 


where 


IR  (13)  is  Az  encoder  angle  found  in  FMRDRM,  Item  13 
Ic  (3)  is  Az  bias  found  in  Calibration  Record  FMAACC,  Item  3 

2.  On  or  After  12  March  1971 


Az  =  IR  (13) 


fThis  R  is  computed  by  the  ARTP  and  only  approximates  the  true  R. 
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G. 


El 


1.  Prior  to  12  March  1971 


El  =  IR  (14)  +  Ic  (5)  -  AE 

where 

IR  (14)  is  El  encoder  angle  found  in  FMRDRM,  Item  14 
1^  (5)  is  El  bias  found  in  FMAACC,  Item  5 
AE  is  tropospheric  refraction  correction 
2.  On  or  After  12  March  1971 


El  =  IR  (14)  -  AE 


H. 


Az  Offset 
Az  offset  (deg)  = 


VHF  Tr  error  (V) 
VHF  LC  sum  (V) 


where 


_ 1_ 

ATr  slope 


1 

*  - 

cos  El 


*  XK 


VHF  Tr  error  is  found  in  FMRDRM,  Item  19,  if  ALTAIR  is  in  the 
point  target  tracking  mode  and  in  FMRDRM,  Item  99,  if  ALTAIR 
is  in  the  chaff  tracking  mode.  ^ 

VHF  LC  sum  is  found  in  FMRDRM,  Item  23,  if  ALTAIR  is  in 
point  target  tracking  mode  and  in  FMRDRM,  Item  104,  if  ALTAIR 
is  in  chaff  tracking  mode. 

ATr  slope  is  found  in  Calibration  Record  FMRR05,  Item  2,  if 
ALTAIR  is  in  point  target  tracking  mode  and  in  Calibration  Record 
FMR5CH,  Items  2-6,  if  ALTAIR  is  in  chaff  tracking  mode. 

XK  is  a  conversion  factor  from  rr.rad  to  deg 


# 


The  mode  is  found  in  FMRDRM,  Item  107. 
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J. 


El  Offset 


El  offset  (deg) 


VHP  El  error  (V)  _ 1 _ 

VHF  LC  sum  (V)  AE1  slope 


where 

VHF  El  error  is  found  in  FMRDRM,  Item  20,  if  ALTAIR  is  in  point 
target  tracking  mode  and  in  FMRDRM,  Item  100,  if  AI  TAIR  is  in 
chaff  tracking  mode. 


AE1  slope  is  found  in  FMRR05,  Item  4,  if  ALTAIR  is  in  point  target 
tracking  mode  and  in  FMR5CH,  Items  12-16,  if  ALTAIR  is  in  chaff 
tracking  mode. 


K.  RCS# 


where 


RCS  =  10  log 


„  (TGC+RGC)/6  " 

2  ) 

B  *  R  *  2 

I 

1 

1 

„  r  2 

j  — 

B_  *  P2 

R 

-5 

) 

(10  *  KSENS) 


B1  is  a  calibration  constant  stored  in  the  program.  For  VHF  RCS, 

Bi  =  0.460  *  10  "13;  for  UHF  RCS,  B1  =  0.21  *  lo"13 

TGC  is  the  TGC  attenuation  (db)  found  in  FMRDRM,  Item  26 

RGC  is  the  RGC  attenuation  (db)  found  in  FMRDRM,  Item  25 

Bj.  is  a  constant  stored  in  the  program  as  a  function  of  waveform. 

It  takes  on  the  following  values: 


Waveform 

Value 

Long  chirp 

10 

Short  chirp 

4.  5 

CW 

1 

CWL 

10 

P  is  peak  power  found  in  FMRDRM,  Item  22  (VHF)  and  FMRDRM, 
Item  32  (UHF).  It  is  used  in  counts. 


U 

The  RCS  computed  by  ALTCEP  should  not  be  used  for  signature  analysis  since 
Bj  is  a  nominal  calibration  constant  and  is  not  measured  each  mission. 
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1^  is  the  frequency  and  polarization  dependent  target  return  found 

in  FMRDRM,  Item  23  (VHF  LC),  FMRDRM,  Item  29  (VHF  RC), 
FMRDRM,  Item  33  (UHF  LC),  and  FMRDRM,  Item  39  (UHF  RC). 
These  items  are  used  in  counts. 

KSENS  is  the  receiver  sensitivity  bit  found  in  FMRDRM,  Item  51. 

it 

L.  VHF  Waveform 

VHF  waveform  is  determined  by  combining  FMRDRM,  Items  3  and  28. 

M.  VHF  Peak  Power 

VHF  peak  power  (MW)  is  found  in  FMRDRM,  Item  22. 

uu 

N.  VHF  Attenuation 

VHF  attenuation  (db)  =  I  (25)  -  (10  *  KSENS) 

K 

where 

I  (25)  is  VHF  RGC  attenuation  found  in  FMRDRM,  Item  25 
K 

KSENS  is  the  receiver  sensitivity  bit  found  in  FMRDRM,  Item  51 

O.  VHF  PRF 

VHF  PRF  =  106/CATF  (5) 

where 

CATF  (5)  is  the  VHF  PRI  (ps)  found  in  FMCATF,  Item  5 

P.  UHF  Waveform^ 

UHF  waveform  is  found  in  FMRDRM,  Item  38 

Q.  UHF  Peak  Power 

UHF  peak  power  (MW)  is  found  in  FMRDRM,  Item  32. 


Called  V  M  in  listing. 

uu 

Called  VHF  GN  in  listing. 
+  Called  U  M  in  listing. 
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*  t 


I 

fl 

1 

#  1 

R.  UHF  Attenuation 

UHF  attenuation  (db)  =  I_  (35)  -  (10  *  KSENS) 

K  i 

l 

where  , 

I  (35)  is  UHF  RGC  attenuation  found  in  FMRDRM,  Item  35 

K  | 

*  KSENS  is  the  receiver  sensitivity  bit  found  in  FMRDRM,  Item  51 

S.  VHF  Range  Track  Mode  j 

VHF  range  track  mode  indication  combines  the  range  track  mode  found 
in  FMRDRM,  Item  62,  and  the  RTC  loss  of  signal  indicator  found  in  FMRDRM,  Item  73. 

> 

The  loss  of  signal,  shown  as  an  L  to  the  right  of  the  track  mode,  indicates  that  the 

signal  is  either  saturated  or  of  insufficient  amplitude.  : 

T.  Angle  Track  Mode 

Angle  track  mode  indication  combines  the  angle  track  mode  found  in 
FMRDRM,  Item  69,  and  the  ATC  loss  of  signal  indicator  found  in  FMRDRM,  Item  74. 

U.  Chaff  Tracking  Gate  Width 

Chaff  tracking  gate  width  is  found  in  FMRDRM,  Item  107:  i 


Code  Width  (km) 


W.  VHF -UHF  Range  Difference 

VHF -UHF  range  difference  is  found  in  FMRDRM,  Item  75. 

X.  UHF  Range  Track  Mode 

UHF  range  track  mode  combines  the  UHF  lock -on  indicator  found  in 
FMRDRM,  Item  125,  and  the  UHF  loss  of  signal  indicator  found  in  FMRDRM,  Item  126. 

Y.  Tag  Time 

Tag  time  (s)  is  found  in  FMCATF,  Items  1  and  2. 


Z. 


UHF  PRF 


UHF  PRF  = 


10 


CATF  (6) 


where 


CATF  (6)  is  the  UHF  PRI  (ps)  found  in  FMCATF,  Item  6. 
AA.  No.  of  ARS  VHF  Targets 

No.  of  ARS  VHF  targets  =  C^^  -  1 


where 


C2^  is  total  no.  of  VHF  targets  being  tracked,  including  VHF  dummy 
target,  found  in  FMCATF,  Item  21. 


BB.  No.  of  ARS  UHF  Targets 

No.  of  ARS  UHF  targets  =  C^  -  C^^  -  1 

where 

C^  is  total  no.  of  targets  being  tra  1  .ncluding  VHF  and  UHF 
dummy  targets,  found  in  FMCATF,  Item  11. 


CC.  ARS  Target  No.  of  ALTAIR  Tracked  Target 

The  procedure  for  estimating  which  ARS  target  is  the  ALTAIR  tracked 
target,  for  each  record  on  the  catalog  tape,  is: 
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1.  Determine  no.  of  FMRDRM  data  blocks  (0,  2,  4,  or  6)  in  the 
record  from  FMRDRD,  Item  2. 

2.  Compute  TAL  and  R  of  ALT AIR  tracked  target  for  each 
FMRDRM  data  block. 

3.  Datermine  no.  of  minor  cycles  (1  to  8)  in  the  record  from 
FMCATF,  Item  10. 

4.  Compute  TAL  for  each  minor  cycle. 

5.  Compare  TAL' s  for  all  minor  cycles  with  TAL’s  for  all 
FMRDRM  data  blocks.  Determine  minor  cycle  and  FMRDRM 
data  block  closest  in  time. 

6.  Obtain  R  for  each  ARS  tracked  target  (identified  in  FMCIDB) 
from  FMCTDB. 

7.  Select  ARS  target  closest  in  R  to  ALTAIR  tracked  target. 


Priority 

Priority  is  found  in  FMCIDB,  Item  4. 


III. 


OPERATION 


A.  Input 

Start  and  stop  times  (GMT) 

Skip  time  (s) 

Scale  parameters  for  plots 

u 

Lift-off  time  (GMT) 

A  sample  input  is  given  in  Appendix  A. 
CARD  1  (18A4) 

(Col.) 

1-72  TITL  Title  for  listing  and  plots 

CARD  2  [2(213,  F7.3),  F7.3,  315,  213,  F7.3] 


1-  3 

ISTART(l) 

(h) 

4-  6 

ISTART(2) 

(min) 

■  Start  time  (GMT) 

7-13 

TSTART 

(s) 

14-16 

ISTOP(l) 

(h) 

17-19 

ISTOP(2) 

(min) 

i  Stop  time  (GMT) 

20-26 

TSTOP 

(s) 

27-33 

DELTAT 

Seconds  between  each  line  output  in  meti'ic  listing 

34-38 

IDELX## 

Ordinate  scale  spacing  (dbsm/in)  (10) 

39-43 

MINXPL## 

Minimum  ordinate  of  plot  (dbsm)  (-50) 

44-48 

IDELT## 

Abscissa  scale  spacing  (s/in)  (5) 

(If  <  0,  no  plots  are  produced) 

49-51 

ILONCH(l) 

(h) 

52-54 

ILONCH(2) 

(min) 

Lift  -off  time 

55-61 

YLONCH 

(s) 

# 


Normally  left  blank  and  lift-off  time  from  Calibration  Record  FMGLOT  is  used. 


If  left  blank,  program  sets  to  indicated  value. 
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B. 


Output 


METRIC  LISTING 
GMT  h,  min,  and  s 
TAL 
Alt 
R 

R  (m/s) 

Az 
El 

Az  offset 
El  offset 
VHF  LC  RCS 
UHF  LC  KCS 

u 

VHF  waveform: 

VHF  power  (dbw) 

VHF  attenuation 
VHF  PRF 

t 

UHF  waveform 

UHF  power  (dbw) 

UHF  attenuation"^ 

VHF  range  track  mode:  T  =  track;  C  =  coast;  L  =  lock-on; 

S  =  slaved 

VHF  range  tracker  loss  of  signal  indication  (L  printed  to  right 

of  track  mode  symbol) 


Called  V  M  in  listing. 

uu 

Called  VHF  GN  in  listing 
^Called  U  M  in  listing. 
^Called  UHF  GN  in  listing. 


CW  =  CW;  S  =  short  chirp;  L  =  long  chirp; 
CL  =  CW  long  pulse 

(db) 
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Angle  track  mode:  T  =  track;  C  =  coast;  M  =  manual;  S  -  slaved 

Angle  tracker  loss  of  signal  indication  (L  printed  to  right  of 

angle  mode  symbol) 

Chaff  tracking  gate  width  (km) 

Track  reference:  TE  =  trailing  edge;  CN  =  centroid; 

LE  =  leading  edge 

VHF  -  UHF  range  difference 

UHF  range  track  mode:  T  =  track;  blank  =  not  in  track 

UHF  range  tracker  loss  of  signal  indication  (L  printed  to 

right  of  track  mode  symbol) 

ARS  TARGET  LISTING 
GMT  h,  min,  and  s 
TAL 

Tag  time  (s) 

VHF  PRF 
UHF  PRF 

No.  of  VHF  targets 
No.  of  UHF  targets 

ARS  target  no. ,  frequency,  and  priority  for  target 

closest  in  range  to  the  ALTAIR  tracked  target 

Absolute  range  difference  between  selected  ARS  target 
and  the  ALTAIR  tracked  target 

ARS  target  test  data  appear  only  when  one  of  the  following  items  change: 
VHF  PRF,  UHF  PRF,  no.  of  VHF  targets,  no.  of  UHF  targets,  ARS  target  no. ,  fre¬ 
quency,  and  priority. 


PLOTS 

RCS  plots  vs  TAL  of  ALTAIR  tracked  target 
(VHF  and  UHF  LC  and  RC) 

Sample  outputs  are  presented  in  Appendix  B. 


12 


IV.  PROGRAM  LIMITATIONS 

Start  time  Must  be  before  end  cf  tape 
Stop  time  Must  be  before  start  of  tape 
DELTAT  Must  be  a  multiple  of  0.025  s 


i 
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V. 


PROGRAMMING 


A.  ALTCEP  (see  Appendices  C  and  D. ) 

ALTCEP  is  the  control  section  of  ALTCEP.  ALTCEP  reads  the  input 
cards,  calls  the  subroutires,  and  prints  the  data. 

B.  CHEAP  (see  Appendices  E  and  F  ) 

CHEAP  is  used  to  process  calibration  and  format  records  which  are 
recorded  before  the  data  records  on  the  ALTAIR  transcription  tape.  CHEAP  lists 
format  and  calibration  records  named  in  the  common  statement.  These  are  stored 
and  unpacked  for  later  use  by  the  main  processing  program.  A  sample  CHEAP  output 
is  given  in  Appendix  G. 

The  minimum  size  needed  for  the  item  array  may  be  calculated  by 
the  following  equation: 

Item  size  =  6  *  (Total  number  of  items  stored)  +  (6  *  130) 

Calling  Sequence:  Call  CHEAP  (*) 

*  =  A  return  point  specified  by  a  statement  number 
in  the  calling  program.  Used  for  aborting  job 
by  main  progra  n  if  wrong  tape  is  mounted. 

CHEAP  calls  the  following  subroutines: 

BREAPS  (entries  BREAPS  and  BREAD);  HPRR  (entries  HPRR 
2 

and  NAMET);  and  FORM.  CHEAP  also  calls  WHICHV,  a  360 -system  subroutine 
indicating  whether  a  job  is  being  run  under  the  time -sharing  (CMS)  or  Batch  (OS)  systems. 

C.  GET  and  IGET 

2 

GET  and  IGET  are  entries  to  subroutine  GETS. 

These  routines  will  locate  any  data  item,  unpack  it,  and  interpret  it 
according  to  the  information  in  the  format  table.  They  will  return  the  item  as  a 
binary  integer  (in  the  case  of  IGET)  or  as  a  floating  point  number  (in  the  case  of  GET). 
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GET  (or  IGET)  requires  three  arguments: 

GET  (Format,  Base,  Item) 

Format  Relevant  format  table  address 

Base  Base  address  of  data  block  desired 

Item  Specific  item  number 

D.  LTIME  (see  Appendix  H. ) 

Subroutine  LTIME  unpacks  the  lift-off  time.  The  call  statement  is 
LTIME  (LOT,  IH,  IM,  IS,  IT). 

INPUT 

LOT  Unpacked  lift-off  time 

OUTPUT 

IH  h 

IM  min 

IS  s 

IT  ms 

E.  GMTTAL,  CATIME,  and  RADART  (see  Appendix  J. ) 

GMTTAL,  CATIME,  and  R/DART  are  entries  to  Subroutine  GMTTAL. 
1.  GMTTAL 

Entry  GMTTAL  converts  start  time  and  stop  time  on  input 
data  cards  to  total  GMT  s  and  computes  TAL.  The  call  statement  is  GMTTAL 
(IGMTH,  IGMTM,  GMTS,  TAL). 

INPUT 

IGMTH  Hours 

IGMTM  Minutes 

Seconds  and  decimal  fractions  of  seconds 


GMTS 


OUTPUT 


TAL  Time  (s)  after  launch 

STORED  IN  COMMON 
TLONCH  Lift-off  time  (total  GMT  s) 

2.  CATIME 

Entry  CATIME  obtains  GMT  from  FMCATF,  unpacks  and 
converts  it  to  total  GMT  s,  and  computes  TAL.  The  call  st  ement  is  CATIME 
(FMCATF,  IADD,  IGMTH,  IGMTM,  GMTS,  TAL). 

INPUT 

FMCATF  FMCATF  format  table  address 

IADD  Base  address  of  FMCATF  data  block 

OUTPUT 

IGMTH  Hours 

IGMTM  Minutes 

GMTS  Seconds  and  decimal  fractions  of  seconds 

TAL  Time  (s)  after  launch 

STORED  IN  COMMON 
TLONCH  Lift-off  time  (total  GMT  s) 

3.  RADART 

Entry  RADART  obtains  GMT  from  FMRDRM,  unpacks  and 
converts  it  to  total  GMT  s,  and  computes  TAL.  The  call  statement  is  RADART 
(FMRDRM,  IRDRM,  IGMTH,  IGMTM,  GMTS,  TAL). 

INPUT 

FMRDRM  FMRDRM  format  table  address 

Base  address  of  FMRDRM  data  blocx 


IRDRM 


OUTPUT 


IGMTH  Hours 

IGMTM  Minutes 

GMTS  Seconds  and  decimal  fractions  of  seconds 

TAL  Time  (s)  after  launcn 

STORED  IN  COMMON 
TLONCH  Lift-off  time  (total  GMT  s) 

F.  REFC  (see  Appendix  K. ) 

The  tropospheric  refraction  correction  subroutine,  REFC,  is  based 

3 

on  tropospheric  refraction  tables  in  PPP-36.  A  modified  version  of  this  subroutine 
is  now  in  use. 

The  call  statement  is  REFC  (E,  R,  DEE,  DRR) 

E  Uncorrected  El  (must  be  between  0°  and  90°) 

R  Uncorrected  R  (ft) 

DEE  El  tropospheric  correction 
DRR  R  tropospheric  correction  (ft) 

The  corrected  values  to  be  computed  after  exiting  from  REFC  are: 

El  =  E  -  DEE 

R  (ft)  =  R  -  DRR 

G.  BZERQ  (see  Appendix  L.  ) 

Subroutine  BZERO  is  necessary  if  a  floating  point  item  is  scaled  BO 

in  a  format  or  calibration  table  description.  BZERO  is  called  after  the  item  has  been 

2 

extracted  by  subroutine  GETS.  BZERO  normalizes  the  data  item  and  puts  the  decimal 

point  in  its  proper  position.  The  call  statement  is  BZERO  (yy). 

INPUT  AND  OUTPUT 
yy  Iter  to  be  processed 
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H. 


I  BIT  (see  Appendix  M.) 


Function  IBIT  identifies  which  bit  is  set  from  a  group  of  bits,  assuming 
one  and  only  one  bit  is  set.  For  example,  track  reference  is  determined  from 
FMRDRM,  Item  61.  There  are  three  bits  of  binary  data,  numbered  from  right  to 
left,  in  the  the  format  ABC  where 

A  indicates  leading  edge 
B  indicates  ^.entroid 
C  indicates  trailing  edge 

The  cal)  statement  is  IBIT  (ITEM,  NBIT). 

INPUT 

ITEM  Designation  of  group  of  bits  to  be  searched 

NBIT  No.  of  bits  in  group 

OUTPUT 

No.  of  the  bit  that  is  set 

J.  METAZ 

METAZ  is  the  plotting  routine. 

K.  Plotting  System  Subroutines 

They  are  PLTND,  STOIDV,  and  REREAD. 
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APPENDIX  A 
ALTCEP  INPUT 


ALTCEP  TEST  RUN 


CARD  1 


?  “  3  ?  s  S  ?  a  ^OMi^nu,0,,0,0  0  0  P 0  °,9  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0 0  0  C  0  0  0  0  0  0  0  0  D  0  0  0  0  0  0  0  n  0  0  n  0  0  0 

I  1  1  1  <  I  1  <  1  ..........  .  ‘”3  4ZS  21  :13  31,  31  32  33  S2SJS<SS  58  5;5B  59  80  6I  62  83  B<C5  66  6/B8£3;0;i  »n;<;W8  ;//3;i80 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1  ill  i  III  I  ill  1 1 1 1 1 1 1 1 1 1 1 1 1 1 : 1 1 1 1 1 1 1 1 1 1 ,,  I ,,,,,  1]M  M  n  ,, , 

2  2  2  2  2  2  2  2  2  22222222222  2  222222222222222  1  2222222222222l2222  2  22222222222222222222222 

3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3 ,3-1  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3 
4444444444444  444444444444444444444444444444444^44444444444444444444444444444444 
5  5  5  5  5  5  15  5  5  5  5  5  5  5  5  5 5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5 


1 11 i 1  7777777777777777777777777777777777777777777777777777777777777 


11111111111111 


8  C  8  3  3  C  0 1 E  S  8  8 


■  "  i  ?  ?  ['  ? 5 11 .  5  5  .^^9999999999999999999999999  J;(99999999999999999?9g'''ig'’q 

’,'tLIS  ?,J2‘  ‘‘  •  30  31  3?  33 11  35  38  52  '.f  39  40IU7O<<<S4M2  <8  49  50  51  j25]S4  51  86  5/  SR  59  63  61  62  6]  (468  (6  62  (869  20  (I  l>n 


6  58  0-300  6  58  30.200  0.1 


CARD  2 


00009000  00  00000000  00  00  00000000  000000000000300000000000000000000000000000 

:  !  1  1  ?  ?  V  ?  ?  V  V  ?  2°  23  22  23  24  25  28  27  28  29  50  31  32  33  34  3S  36  J7  39  3,  4° 41  42  <3  88  25  <8  42  89  6,  50  51  53  54  55  56  $7  SI  59  SO  61  62  63  68  68  88  67  68  69  70  /I  82  73  74  25  26  2/  28  '9  I] 

1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  I  I  1  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  1  I  |  |  |  |  |  |  |  1  |  1  |  |  |  |  1  |  1  |  1  |  |  |  |  |  |  |  |  |  |  |  ]]  |  ,,,  , 

22222222222222222222222  22222222222222222222222222222122222222222222222222222222 
333313333  333333333  3  333333  33333333333333333333333333333333333333333333333333 

44444444444444444444444444444444444444444444444444444444444444444444444444444444 
5555  355555555555  55555555555555555555555555555555555555555555555555555555555555 
0  5  0  0  6  0  6  6  6  0  6  0  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  S  6  6  6  6  6  6  6  6  6  6  6  6  6  6  6  0 6 

7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  M  7  7  7  7  7  7  7  7  7  H  /  7  /  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  /  7  7  7  7  J  7  7  7  7  7  7  7  7  7  7  7  7  7  /  7  7  7  7 

89836  888  88888808  808  888080  88808888888888888888088088888888888888888868086800 
9  5 1 5  9  9  9  8  9  9  0  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  T  '  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  9  n  g  Q  g  •,  q  n  Q 

■■  k  '  3  ^  15  IS  17  is  13  2021  22  23  24  26  <5  27  28  29  33  <•  3;  i2  "<  35  36  32  'JB  39  40  41  42  63  44  45  46  62  33  43  5Q  M  52  S3  S4  55  56  22  58  S3  60  6!  52  63  61  G5  66  62  68  66  73  21  22  7*  .4  75  4  V  23 1-  S3 
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APPENDIX  B 
ALTCEP  OUTPUTS 
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APPENDIX  C 

ALTCEP  PROGRAM  LISTING 


COMMON/ 8EAD/LN, IFLG, I  ADD , FHRD l D , FMCATF , FMCS AD , FMCMD8 » FMCT I  8 , FMC I DB 

1, FMCTDB,FMRDRD,FMRDRM,FMRDRT,FMGL0T,FMRR05,FMAXSP,FMBIAS,FMR5CH 

2 ,  r  MRCHF ,FMAACC,FMRR11,NAME( 19),NII18),IX(18),1 TEM( 8000) 
COMMON/TITLE/ IPARRI 10 ) 

COMMON  /LAUNCH/  TLONCH 

INTEGER  ♦  2  I  TP M , TKREF1 , TKREF , L8 I T , I TL B I T , I RTB I T 

INTEGER  ♦  2  I25MIL,  IAM8,  I8UFF,  I8UF0,  IFLAG,  IGMTH,  IGMTH2 , 

1  IGMTM,  IGMTM2,  ILONCH,  IPAGE,  I  PR  I ,  IPRI02,  (START,  ISTCP,  ITAR2, 

2  ITARG,  ITKMOD,  ITYPE2,  IUPRF2,  IVPRF2,  JTYPE,  KFLAG,  KCN1C,  LCS, 

2  LOSANG,  LOSRAN,  MAX8UF, IGATE  « 

3MAXL 1 2 ,  MAXLIN,  MFLAG,  MINOR,  MODE,  M0CTK1,  M0DTK2 ,  N25MIL, 
ANCHECK,  NFL  AG ,  NRDRM,,  NTGRN,  NTGRO ,  NTCT ,  NTOTE ,  NTYPE,  NUHF2, 
5NUHF0,  NUMBUF,  NUMLIN,  NUMREC,  NUMUHF ,  NUMVHF ,  NVHF2,  NVHFC, 
6UGAIN,  UMOOE ,  VGA  IN,  VMODE 

INTEGER  *  2  N YR ( 2  I , M TH 1 ( 2  I , MTH 2 ( 2 ) , MTH3 ( 2 ) , MTHA ( 2 ) , NDAY 1 ( 2 ) 
1*MTH(2)*IYR(2)«MTHA(2)*NDAY(2)«LA8YAX(9»A) 

REAL  ♦  8  GMT S ,  GMTS2,  TAGTI2,  TAL,  TAL2,  TIME,  TMINOR,  TRTAL2, 

IT  S  TAR  T ,  T S TOP 
REAL  ♦  8  TLONCH 

DIMENSION  CROSFCI 1210, A),  DRANG2I570),  GATSEL ( 6  )  ,  GMTS2I570), 

1  I AMB ( 20 ) , L  B I T ( 3  1, 

I  I  AX ( 2 ) ,  I8UF1C20A8),  I BUF2 ( 204  8  )  ,  1 8UFF ( 5 70 , 7 ) ,  I8UF0I7), 

2 1 GMTH2  (5701,  10MTM2<570),  IL0NCHI2),  INSTRI16),  IPRK20), 

3  I  PR  I  02 ( 570  J ,  I S  T  ART ( 2 ) ,  I S  TOP ( 2 ) ,  ITAR2I570  I, 
3ITARG(20),IXTICK(9),KXTICK(7),JXTICK(9),ITICK(9),ILCT(3), 

AIT  I  ME ( 1210 J ,  I TKMOD (21,  ITYPE2I570),  IUPRF2I570),  IVPRF2I570), 

5  I  XPOS  (  2  J ,  IYPnsm,  JAX(2),  JCROSSI  1210,  A)  ,  JTICKI9),  JTYPE  (  3  ) , 
6KAX(2),KTICK(9) , L  A8EL ( 16) ,LAX (21  , LOCX ( 16) , 

7  L0CYU6I,  LOS ! 2 ) *  MES  AGE  (18),  M0DEI5),  MOD  TK  1  (  A  )  ,  MCDTK2  (  A  )  , 

7  NTYPE ( 256 1 , 

8NUHF2I 570) ,  NVHF2I570),  PLOT  I M (  1 2) 0  I ,  RANGE  2 (10),  TAGTI2I570), 

9  TAL2I  570  ),  TRTAL2I  10),  XLEFTI12I0),  XRITEI1210),  XXMES(12), 

A  Yl.EFTI  1210),  YRITE11210) 

DIMENSION  MC A TF ( 8  1  ,  TKREF 1(3) 

EQUIVALENCE  (IBUFFI1,)),  IVPRF2I1)),  I  I  BUFF  11,2),  IUPRF2 ( I  )  )  , 

1  (  I8UFFI  1,31  ,  NVHF  2  (  1 )  '/ ,  (IBUFF(1,A),  NUHF2I  1 )  )  ,  «IBUFF(1,5), 
2ITAR2II)),  (  I  8UFF  (  1 , 6  I  ,  I  TYPE  2  «  l »  )  ,  I  I  BUFF  11,7),  IPRI02UM 
EQUIVALENCE  (MINXPL,  LABEL  (  8  )  )  ,  (PLOTIMU),  ITIME(D), 

1  ( CROSECI 1, 1  I ,  JCROSS (1,1),  XLEFT(l)),  (CR0SECI1,2),  XRITE(l)), 

2  (CROSECI  1,3)  ,  YLEFT(U),  ( CROSEC  (  1 ,  A  ) ,  YRITE1D) 

EQUIVALENCE  ( FMCATF  ,  MCAYF (1)1 

EQUIVALENCE  ( MVR , NYR ( 1  I) , 1  NTH , MTH( 1 ) ) , I NTHA , MTH A ( 1 ) ) , 

1 ( MDAY , NDAY ( 1 ) ) 

DATA  lYR/«  71*/,MTHl/»  JA’/,MTH2/»  FEV, 

IM/ Z5FFFFFFF/ , M  TH  3/ 1  MA'/,MTHA/»  OR  * / , KK / 1 677 72 1 6/ 

2.L/Z7FFFFFFF/ , NDAY  l/*  1 2 */, MM/ ZOOOOF 080/ , NN/ 1 00008000/ 

DATA  CDN10  / 10./,  C0N30  /30./,  CONALT  /.30A80061/,  CCNDEG  /360./, 
1C0N0FP  /. 078125/,  CONPDW  /l.E-10/,  CONRAN  /  .  00091 AAOl  8 3/ , 

2C0NRAT  /. 0027905329/,  CONRCS  / . 55  1  1 8 1 1 C2A/ ,  CONUP  /2.666666G7/ , 

3  CLNVP  /l. 333333333/,  GATSEL  /O. ,  2.5,  5.,  10.,  20.,  AC./, 

3  I2T22  /A19A30A/, I AX/102C, 1020/, IFL/O/, 
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5ITICK  /  9*  4C0/,  IXPOS/  1911,770  /,IYPCS/  175.1309/ 

6JAX/400,400/,JTICK/9*325/,JTYPE/’  ’,’  V’,’  UV. 

8 1020 1330 l°?#>4n  To‘°  /10/’  K°NTAG  /lC00000'’  “TICK  MOO,  710, 
81020,  1330,  1640,  1950,  2260,  2570,  2880/,  LABXAX  /'TAL  •/ 

oata  l  a  b  ya  x  /  *  vh  • ,  •  f  -,’lc’,’  r.,.Cs’,’  (•  •obs*sm»,m  ».'vh* 

( * , 'ob* , •  sm • , * ;  vu-.-  rv^*, 

2  *  ’  DB  ***UH*,'F  •  »  *RC  •  ,  •  R’,*CS*,’  ( • , ’OB’ , ’SM’ , • )  ’/ 

IA7l;5//^c°rvfl/7/IL^cx/963, 1273,  1 583f  1893,2203,25 13’2823* 

2  »7n  !.25?’  ??°!  7I°’  1020’  l330’  ‘MO.  1950,  2260, 


2  2570,  2880/,  LOS  /’  •,  *L  »/,  MAX  /8192/,  MAXBUF  /570/, 


*  -  *  w  »  f  HMA  /Oi7£f  , 

3  MAXLIN  / 55/ ,  MFLAG  /!/,  *QDE  /’CW,  •  S*, 


4  MOOTK 1  / • T  • ,  *c  • ,  'L 


'CL 


C 

6000 

6001 

6010 


—  /, 

^NFLAG^I/^NTOrS  /O/^  “  X  U  HARD  CCPY  LENS'AN0  CU  T .  "iH  ANKS  •  / 

8JXt!ck/9*945/33°’ l6A0’  1950 *2260* 25 70, 2880/, I  XT ICK/9* 1020/, 

DATA  ICMS  /'CMS  * / , LB  I T /  •  »,'T  ',*L  */ 


6020 

6030 


6040 

6050 

6060 

6070 

6080 

6090 

6100 

6110 


FORMAT  (2(213,  F7.3),  F7.3,  3 1  5 , 2 1  3 , F 7 . 3 ) 
FORMAT(2(2I3,F7.3),F7.3,213,F7.3> 

,F°R?AT  |6?**G?!!’8X,’TAL'  ’5X,’  ALT',  2X,'  RANGE  ',3X,*  RANGE  ’,4X,’A2*, 
13X,  FL  *3X,  M  , AX , *  EL • , 4X , * VHF • , 4 X , • OFF ' , 2X ’ V ’ , 2X • • VHF •  IX 

i;rp;RiMr^;fx;,u,’2x,uHF,’ix’,uHF,’2x*,MooE;'2x',cHFj, 

?JX*  ;.lX;  VR’,3X,'UHF'/38X,  'RATE*,  15X,  'OFF',  3X,’CFF',3X, 

4  LC  ,5X,  LC  , 3X, ’M’ , 2X , ’PWR • ,2X, *GN’ , 1 X, * PRF’ ,2X, • M * , 2X. • PRR* 

52X..GN'.8X,-GrF.,lX,'REF.,2X,'-UR.,i;,.RTK'/i;^.(SE’f^;3^  ’ 

6  (KM)  ,2X, ’ (KM) • , 3X, ’ ( M/S) ’ ,4X, • <D) ’ ,2X, • (D) ’ ,2X, ’ ( C) • 

I  »X »  !nli|  *  *  lX*  ’ ^  OBSM )  '  ,2X,  '  (  CBK)  •  ,  IX ,  ’(DB  )  '  ,6X, 

,™  /  X'  ,Dn*  •  1X’  R’  »2X,  ’A  ’,  IX,  ’  (  KM)  •  ,5X,  ’  (KM)  •  /  ) 

1U?I5  A2  F5  I‘?IF?ni25I;f,«*il,i6*l,F5,l,2F6‘2,2F6'1,lX**2»F5-1* 

1cnDu»TAf,F5*1,IA,2lX,2AlE5*1’lX’A2’lx’F5*2*2X,2Al) 

F°RMAT.  L  .  6X’  ,GMT’'  l0X’  5  X  ,  ’TAG  TIME’,  5X,  'V  PRF’, 

218^  V,«r  J  I:  '??*  9X’  POSSIBLE  SYLVANIA  TARGET*  / 

218X,  (SEC)  ,  6X ,  ’(SEC)’,  24X ,  ’VHF’,  3X,  ’UHF ’ ,  4X,  ’ARS  TARG  • 
33X,  ’PRIOR.’,  3X,  ’RANGE  OIFF.  (KM)»  /)  *  ’ 

FORMAT  I//3’  Fuo3infF1.1:3’  21 9’  I7f  ,6’  I9’  1X’  A2  ’  ,8’  FM-2) 

FORMAT  (//•  UP  TO',  15,  •  RECOROS  SKIPPED’//) 

FORMAT  (1H1) 

FORMAT  (’  START’, 213, F7. 3,  •  STOP • , 2  I  3 , F 7. 3 , •  DELTA’, F7. 3. 
1F0RMAt’(16A4)1ELTA  X,’I5,,  °ELTA  T’’I5»'  LAUNCH’,  213,  F7. 3) 

FORMAT  (1614) 

FORMAT  (18A4) 

FORMAf  (  IH 1 ,  1«A4,  48X,  ’PAGE’,  15) 


IGLOT(K)  =  I G  E  T ( FMGLOT , I ADD,K ) 
RRll(K)  =  GETIFMRRU,  IADO,K) 

AACC(K)  =  GETIFMAACC, I  ADO , K ) 
CATOAT(K)  =  GETlFMCATr, I ADO,K / 
ICAOAT(K)  =  IGLTIFMCATF, IAOO.K) 

I  FOR 0 ( K  )  =  IGETIFMRORO , I RORD,K) 
IRCOAT(K)  =  IGFT(FMRDRM,IROPM,K) 
ROOAT(K)  =  GETIFMRORM, IRORM,K) 
lOBCAT(K)  =  IGFTIFMCIOB,  IOBASE,  K) 
C2(K)  =  GET ( FMR5CH, I ADO,K ) 
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CALL  WHICHVI ID ) 

WRITEI6.1)  ID 

I  FORMAT ( •  I D=  *  A4) 

I  OUT  =  6 

IFI ID.EQ. ICMS) 10UT=8 
C 

XK=. 05729578 
fUUHF  =  0.2099581E-13 
B1VHF  =  0.4604345E-13 
KFLAG  =  l 
IPAGE  =  0 
CALL  CHEAD  IC147J 
C 

READ  (5,  6100)  MESAGE 

READ! 5,6000)  I  5  TART ( 1) , I  START ( 2 ) , T STAR T , I  STOP  I  1 ) , I  STOP (2) , TSTOP, 
I  DEL  TAT, IDELX.MINXPL, IDELT, ILONCHIl) , ILCNCHI 2) , YLONCH 
IF ( ( ILCNCHI I ).NE.O).OR. { I LONC H I  2 ) . NE . 0 ) . OR.  < YLONCH. NE . 0. ) ) 

1  GO  TO  21 
□  0  1 1  I  =  1,3 
ILOTI I )  *  IGLOTI I ) 

II  CONTINUE 
TLCNCH  =  0.0 

IFI ILOTI 1)  .NE.  2)  GO  TO  22 
CALL  LTIME  (ILPT,  NH,  NM r  NS,  NT) 

ILONCHI  l)  *  NH 
ILCNCHI2)  =  NM 

YLONCH  =  NS  ♦  FLOATINT)  *  .001 

TLCNCH  =  ILONCHIl)  *  3600  ♦  IL0NCHI2)  *  60  ♦  YLONCH 
PRINT  23, TLONCH, ILONCHIl), IL0NCH12)  ,NS, NT 
23  FORMAT  I • l'F9.3,3X, 21 12, *0* ), 12, *. » 13, ’  LAUNCH  *  ) 

GO  TO  22 

21  TLCNCH=ILONCH< 1 )*36 00 ♦ILONCHI 2)*60* YLONCH 
22  CCNTINUE 
C 

N  =  IDARRUO) 

J  =  IDARRI9) 

J1  =  MASKIJ,!.) 

N 1  =  MASK  I N, M ) 

J3  =  I  CARR  I  7 ) 

N  3  =  IDARRI8) 

IPROC  =  0 

MYR  =  M0DIN1/256, 65536) 

NTH  =  MODI  J1 ,6*^536  ) 

N2  =  MCDIN1/KK, 256) 

NTHA  =  MAKE  1 1 N ? , MM  J 

MDAY  =  MODI Jl/65536, 65536) 

MD AY  =  MAKE 1 1  MPA Y , NN ) 

IFINYR(2).LT.IYR(2) )G0  TO  12 
IFINYR!2).GT.IYR(2)  )  GO  TO  13 

IF(MTH(2).EQ.MTHl( 21.0R.MTHI2).EQ.MTH2<2) )GC  TO  12 
IFIMTHI 2).NE.MTH3( 2) )G0  TO  13 
IF(MTHA(2).NE.MTH*,<2)  )G0  TO  13 
IFINDAY(2).LT.N'JAYl(2i  )  IPROC  *  1 
GO  TO  13 

12  IPROC  =  1 

13  CCNTINUE 

PRINT  991, IPROC, J,N,  3,N3 
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rv’ 


n  <e  d* 


o  o  o 


991  FORMAT  1  /  /  *  IPRPO',12//*  ARTP  VERSION  CATE  *  *,2A4/, 

1»  RTS  VERSION  PATE  =  *,2A4//| 

INITIALIZE  PLOTTING 

NPLCT  =  0 
PLMINX  =  -50. 

IF  ( IDELX )  6,  %  2 

2  PLMINX  =  MI  NX PL 
GO  TO  6 

3  M I  NX  PL  =  -50 
IDELX  =  10 
IDELT  =  5 

6  WRITE  (6,6070)  I  START  Ill, 

1 1  START ( 2  )  »  TSTART,  ISTCP(l),  IST0PI2),  TSTOP,  DELTAT,  MINXPL, 

2 IDELX,  IDELT, ILONCHI 1 ) , ILONCHI 2 ) , YLONCH 
C 

IF  ( IDELX  .LT.  0)  GO  TO  8 
DELSIX  =  6.  ♦  IDELT 
CALL  STOIDV(MESAfE,71,OI 
CALL  REREAD  (99) 

CALL  METAZ  (20,  0,  0,  0,  0,  0) 

CALL  METAZ  (50,24,2,0,0,0) 

CALL  METAZ  (51,35,0,0,0,0) 

CALL  METAZ  (51,  36,  0,  0,  0,  0) 

CALL  METAZ  (51,  37,  4095,  0,  0,  0) 

CALL  METAZ  (51,  38,  3071,  0,  0,  0) 

CALL  METAZ  (51,  31,  0,  0,  0,  0) 

CALL  METAZ  (51,  32,  512,  0,  0,  0) 

CALL  METAZ  (51,  33,  4095,  0,  0,  0) 

CALL  METAZI51,  34,  3583,  0,  0,  0) 

CALL  METAZ  (50,  5,  2,  0,  0,  0) 

CALL  METAZ  (51,  27,  39,  0,  0,  0) 

CALL  METAZ  (50,  3,  l,  0,  0,  0) 

CALL  METAZ  (5,  0,  720,  1600,  72,  MESAGE(ll) 

WRITE  (99,6001)  I  ST  ART (  1 ) , I  ST AR T ( 2 ) , T S T AR T , I  STOP ( 1 ) , I S TOP ( 2 ) , 

ITS  TCP, DELTAT, ILONCHI 1 ) , I LCNCH ( 2 > , YLONCH 
READ  (99,  60801  XXMES 

CALL  METAZ  (5,  0,  720,  1400,  48,  XXNESIl)) 

CALL  METAZ  (13,  l,  0,  0,  0,  0) 

8  CONTINUE 
C 

IF ( ( I  START! 1 ).FQ.O).AND.( I  START ( 2 ) . EQ. C ) . AND. (TSTART 
l.EC.  0.)  )  GO  TO  10 

CALL  GMTTAL  (ISTART(l),  I  START ( 2 ) ,  TSTART,  TSTART) 

C 

10  IF  (  (  I S  TOP (  1  )  .NE.  0)  .OR.  ( I  STOP ( 2 1  .NE.  0)  .OR.  (TSTCP  .NE.  0) 
1)  GO  TO  20 
TSTOP  =  1000000. 

GO  TO  25 

20  CALL  GMTTAL  (ISTOP(l),  IST0PI2),  TSTOP,  TSTOP) 

C 

25  TIME  =  TSTART 
LN  =  l 

NUMLIN  =  MAXLIN 
MAXLI2  =■  MAXBUF  /  10 


C  DETERMINE  TICK  LABFLS 
C 

IF  ( I  DEL  X  .LT.  0)  GO  TO  32 
J  *  TSTART  /  IDELT 
LABEL (  II  =  J  ♦  IDELT 
00  2B  I  =  2,  7 

LABEL (  I  I  =  L ABFL ( I  - L }  ♦  IDELT 
2B  CONTINUE 

TIMMIN  =  LABEL ( 1 1 
TIMMAX  =  LABEL ( 7  I 
00  30  I  =  9,  16 
LABEL (  I  I  =  L ABFL ( I - 1 1  ♦  IDELX 
30  CONTINUE 

PLMAXX  =  LABEL ( 16) 

RCSOIF  =  PLMAXX  -  PLMINX 
WRITE  (99,  6090)  LABEL 
READ  (99,  60B0I  LABEL 
32  CONTINUE 

IF  ( MCATF (  i  !  .CQ.  0)  CALL  EXIT 
00  33  I  *  5,8 

IF  ( MCATF (  I  )  .FQ.  01  CALL  EXIT 
33  CONTINUE 
C 

35  CALL  BREADS  (LN,  IBUF1,  IBUF2,  MAX,  IFL,  INDX,  LEN,  IFLG, 
NUMREC  *  0 
40  CALL  BREAD  (LN) 

GO  TO  (60,  45,  50,  50),  IFLG 
45  NUMREC  =  NUMREC  ♦  1 
GO  TO  40 
50  NFL  AG  =  2 
MFLAG  =  2 
NUMBER  =  NPLQT 

IF  (NUMREC  .EQ.  Cl  GO  TO  142 
IF(NJML IN.GT.50)  WRITE! I  OUT, 6060) 

WRITE! I OU T , 6050 )  NUMREC 
GO  TO  142 
C 

60  IRORO  =  I C  AO A  T ( 191  *  3  ♦  IAOD 
IR0RD2  =  IRORO 
IREL  =  (  I FORC ( 11-11*3 
N25MIL  =  I FORO <  2 ) 

I25MIL  =  I FORO ( 3  I  *  3 
IREL  =  0 
N25MIL  =  2 
I25MIL  =  120 
C 

CALL  CATIME  (FFCATF,  I  ADD,  I GMTH2 ( NUMBLF ) ,  IGMTM2 ( NUMBUF ) 
1GM  TS2 ( NUMBUF I ,  TAL 2 ( NUMBU F  )  ) 

I VPRF2 ( NUMBUF I  =  KONTAG  /  ICADAT(5) 

I UPRF2 ( NUMBUF I  =  KONTAG  /  ICADAT(6> 

C 

GO  TO  (70,1301,  NFL AG 
70  IRORM  =  IRORO  ♦  IREL 

IF  (IRORM  .EC.  IREL I  GO  TO  130 
I50MIL  =  I C AOAT (13)  *  3 
NRORM  =  I C AO A T ( 12  1  *  N25MIL 
K25MIL  =  0 
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DO  120  1*1,  NRCRM 

- 1,11  -ie-  *•* «°  'o  „ 

GO  TO  9  0 

”  - « „ 

TIHE  =  TIME  ♦  OEL  TAT 

GO  TO  72 
CONTINUE 

JfUGA1  2LE*  TSTGPI  G0  TC  75 
number  *  NPLOT 

IF  (NUMREC  .EO.  0)  GO  TO  130 

Jr  ?Mrnn;GI;  =  °’  WR  1  *E  «  IOUT  ,6060) 

»R  I  TE ( IOUT < 6050 )  NUMREC 

GO  TO  130 


LBITI JTLBIT 


”mf :  \ziv,  rMo“3,,m°- 

IRTBIT  *  LBITI IRT6IT) 

JTLBIT  *  IRC0ATI126) 

ITLBIT  =  LBITIJTLBIT  ♦  2) 

IF( IROOAT ( 126). EC. 0) I  TIBI T  *  LBITIJTLBIT  ♦  it 

IFIIPR0C.EQ.1I  GO  TO  711  LHI"JTLBIT  ♦  1> 

PA^at  B  I  R  OOA  T  (  85  )  ♦  RD0AT(86!)  ♦  3.D 
RANRAT  =  ROOA  T(  89  )  *  CONRAT 
EL  *  RDDAT(19)  *  CONOEG 

c*o  .'o”?r'13’  *  cdnoec 

AZ  =  ROOAT (  13)*CCN0EG-»AACC(3)*57.295 
,r  oRG°AT(^»*C0NDEGAAACCr5»*57:295 

RANrFREF?plMr’  RANGE*  OELEL,  DELRNG ) 

CL  -  EL  -  DElIi ~  0ELRNG,  *  -00030A8 
EL  =  EL*. 0179533 
ERAO  *  6378.195 

itf :  H;:««r’e“2*E"'0**2*2‘"‘NGE*E"jc*siN|EEn 

EL  =  EL/. 0179533 
I GATE* IROOAT (107) 

GATE  =  GAT SEL (  IGATE  ♦  I ) 

B2-  IROOA  T ( 109) 

XB2= IROOAT ( 23 ) 

IF(XB2  .N J.  0)  GO  TU  77 
AZOFF  =  10. 

ELOFF  =  10. 

GO  TO  78 

I F ( IGATE  .NE.  0)  GO  10  771 
X9  =  GE T (FMRR05 .  I AOO, ^ ) 

X2  =  GET( FMRR05, IADO, 2 ) 

CALL  BZCRD(X9) 

CALL  B  Z  L  R  0  ( X 2 ) 

I  ROD A  T ( 19)/(XB2*C0S(EL*. 0179533) ) )  *  ( XK/X2) 
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EL  OFF  = ( IROOAT ( 20) /XB2 ) * ( XK/X4 ) 

GO  TO  772 

IF  ( B2  .NE.  0)  GO  TO  773 
AZOFF  =  10. 

ELOFF  =  10. 

GO  TO  78 
773  XC2=C2 (  IGATE-s-1) 

CALL  8ZER0 ( XC  2  I  ' 

AZOFF=  (  IR00AT(99)/(82*C0S(EL*. 0174533) ) ) ♦ ( XK/XC2 ) 
XC4*C2( I  GATE ♦ 1 1  ) 

CALL  BZER0(XC4) 

ELOFF=< IRDDA7(100)/B2>*(XK/XC4) 

772  CCNTINUE 

78  IF  ( I ROOAT ( 3 )  .EQ.  0)  GO  TO  781 

VMOOE  -  MCCEI5) 

GO  TO  782 

781  VMOOE  *  IRDOAT ( 28  ) 

VMOOE  =  MOOEI I8IT(VM00E»4)  ) 

782  VPOWER  =  ROOAT ( 22 )  *  CONVP  *  CONOFP 
IF  (VPOMER  .LE.  C.)  VPOWFR  *  CONPOK 
VPOWER  *  CON  10  *  ALOGIOI VPOWER)  ♦  60. 

KSENS  =  I ROOAT (51)  *  KONIO 
J  =  I ROOAT ( 25 ) 

VGAIN  *  C  J/2  ♦  I J  -  ( J/2 )  *  2)  )  *  3  -  KSENS 
iMOOE  =  IROOAT ( 38  ) 

Ji’OOE  =  MODE (  I  0  l  T  I  UMOOE  1 4 )  ) 

UPOWER  =  ROOAT  (32)  *  CONUP  *  CONOFP 
IF  (UPCWER  .LE.  0.)  UPOWER  «  CONPOW 
UPOWER  *  CONIC  *  ALOG 10 1 UPOWER )  ♦  60. 

J  *  IRODATI35) 

UGAIN  =  I  J/2  ♦  (J  -  (J/2)  *2)  )  *  3  -  KSENS 
KSENS  =  -  KSENS  ♦  10 
TKREF  =  I ROOAT (61) 

TKREF  *  TKREFlt 1 8  I T ( TKREF , 3 ) ) 

ITKMCO(l)  =  IROOAT ( 62  ) 

ITKMOD(I)  *  M00TK1I ITKMOOI  l)  ) 

LOSRAN  =  I ROOAT (73)  ♦  I 
LOSRAN  *  LOS( LOSRAN ) 

I TKMOO ( 2 )  *  IROOAT ( 69 ) 

I TKMOO ( 2 )  *  M0OTK2( IBIT( ITKM00I2) ,4) ) 

LOSANG  *  IROOAT (74)  ♦  I 
LOSANG  *  LOS ( LOSANG ) 

83  =  IROOAT ( 85) 

84  =  IROOAT ( 29 ) 

J  =  IROOAT ( 28) 

B5  =  1. 

IF(J.EC.O)  85  =  10. 

IF  (J  .EQ.  2)  85  =  4.5 
IF  (J  .EQ.  4 )  85  *  10. 

Y  *  IROOAT( 22) 

IF  (  Y  .GE.  8.)  GO  TO  80 
VLEFT  =  PLMINX 
VRITE  =  PLMINX 
GO  TO  81 
80  Y  =  SQRT(Y) 

RGC  =0. 

J  =  1RD0AT(25) 


r>  o  o 


IF  (U  -  (J/2)  *  2 )  .GT.  0>  RGC  *  3. 

RGC  *  RGC  ♦  I  J/2 )  *  3 
TGC  =  I RDDAT ( 26 )  *  3 

VRITE  =  (B1VHF  *  IB3  **  2)  *  (2  **  ( ( TGC+RGC ) /6) )  /  (B5  *  Y))  **  2 
VLEFT  =  VRITE  *  (  XB2  **2) 

IF  (VLEFT  .GT.  0.1  GG  TO  803 
VLEFT  =  PLMINX 
GO  TO  805 

803  VLEFT  =  10.  *  ALOG 10 ( VLEFT I  ♦  KSENS 
805  VRITE  =  VRITE  ♦  (BA  **  2) 

IF  (VRITE  .GT.  0.)  GO  TO  807 
VRITE  =  PLMINX 
GO  TO  808 

307  VRITE  =  10.  *  ALOGIOI VRITE)  ♦  KSENS 

808  B2  *  I RDDAT ( 33 ) 

84  =  IRDDATI 39) 

J  =  I RDOAT I  38  ) 

B5  =  1 . 

IF  (J  .EQ.  2)  85  *  5.01 
IF  (J  .EQ.  A)  85  -  11.21 
Y  =  IRDDAT ( 32 ) 

IF  (Y  .GE.  8.)  GO  TO  809 
ULEFT  =  PLMINX 
URITE  =  PLMINX 
GO  TO  81 

809  Y  =  SORT ( Y  ) 

RGC  =  0. 

J  =  IRDDATI 35) 

IF  ( ( J  -  (J/2)  *2)  .GT.  0)  RGC  -  3. 

RGC  =  RGC  ♦  (J/2)  ♦  3 
TGC  *  IRDDATI  36)  *  3 

URITE  =  (B1UHF  *  (B3  *♦  2)  ♦  (2  **  ( ( TGC^RGC ) / 6 ) )  /  (B5  *  Y ) )  **  2 
ULEFT  =  URITE  *  (  82  **  2) 

IF  (ULEFT  .GT.  0. )  GO  TO  810 
ULEFT  *  PLMIN.'; 

GO  TO  811 

810  ULEFT  =  10.  ♦  ALOGIO ( ULEF  T )  +  KSENS 

811  URIIE  =  URITE  *  (BA  **  ?) 

IF  (URITE  .GT.  0. )  GO  TO  812 
URITE  -  PLMINX 

GO  TO  81 

812  URITE  =  10.  *  ALOGIO(URITE)  ♦  KSENS 
SAVE  DATA  TO  BE  PLOTTED 

81  IF  ( IDELX  .LT.  0)  GO  TO  9A 
NPLOT  =  NPLOT  +  1 

PLOT  I M ( NPLOT )  =  TAL 
IF  (VLEFT  .LT.  PLMAXX)  GO  TO  82 
XLEF  T ( NPLOT )  =  PLMAXX 
GO  TO  86 

82  IF  (VLEFT  .GE.  PLMINX)  GO  TO  8A 
XLEFT ( NPLOT )  =  PLMINX 

GO  TO  86 

8  A  XLEFT ( NPLOT )  =  VLEFT 
86  IF  (VRITE  .LT.  PLMAXX)  GO  TO  88 
XR  I  TE ( NPLOT )  =  PLMAXX 
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GO  TO  881 

88  IF  (VRITE  .GE.  PLMINXI  GO  TO  90 
XRITE(NPLOT)  *  PLMINX 
GO  TO  881 

90  XR  ITE ( NPLOT I  *  VRITE 

881  IF  (ULEFT  .LT.  PLMAXX)  GO  TO  882 
YLEFT ( NPLOT I  =  PLMAXX 

GO  TO  88 A 

882  IF  (ULEFT  .GE.  PLMINX)  GO  TO  883 
YLEFT ( NPLOT )  =  PLMINX  . 

GO  TO  88 A 

883  YLEFT ( NPLOT )  =  ULEFT 

88A  IF  (URITE  .LT.  PLMAXX)  GO  TO  885 
YR I TE ( NPLOT I  *  PLMAXX 
GO  TO  92 

885  IF  (URITE  .GE.  PLMINX)  GO  TO  886 
YRITE(NPLOT)  *  PLMINX 

GO  TO  92 

886  YRITE(NPLOT)  =  URITE 

92  IF  (TAL  .LE.  TIMMAX)  GO  TO  9A 
NUMBER  =  NPLOT  -  1 
IF  (NUMBER  .EQ.  0)  GO  TO  1AA5 
GO  TO  1A? 

9A  TIME  =  TAL  ♦  OTLTAT  -  .001 

IF  (NUMLIN  .LT.  MAXLIN)  GO  TO  96 

95  IPAGE  =  IPAGE  ♦  1 

WRITE( I0UT.611C)  MESAGE, IPAGE 
WRITE ( I UUT  « 601 0 ) 

NUMLIN  =  1 
GO  TO  100 

96  NUMLIN  =  NUMLIN  ♦  1 

100  IF  (NUMREC  .EQ.  0)  GO  TO  105 
NUMLIN  =  NUMLIN  ♦  5 
IF  (NUMLIN  .GT.  MAXLIN)  GO  TO  95 
WRITE( IUUT.6050)  NUMREC 
NUMREC  *  0 

105  W1'  (I0UT,6020)  IGMTH, IGMTM.GMTS, TAL, ALT, RANGE, RANRAT.AZ, 

)  '  ZOFF,  ELTFF ,  VLEFT ,  ULEFT,  VMOOE,  VPOWER,  VGA  IN, 

2IVPw  'NUMBUF ) ,UMOOE,UPOWER,UGAIN, ITKMCDI 1)  , 

3  LOSRAN,  I TKMOH ( 2  )  ,  LOSANG,  GATE  , TXR E F , OEL TR A ,  I RTB I T , I TL 8  I T 
110  K25MIL  =  K25MIL  +  l 

IF  (K25MIL  .EQ.  N25MIL)  GO  TO  115 
IRDRM  =  IRORM  +  I25MIL 
GO  TO  120 
115  K25MIL  =  0 

IRCRO  =  IRORO  +  I50MIL 
IRORM  =  IRORO  +  IREL 
120  CONTINUE 
C 

130  IF  (MFLAG  .EQ.  1)  GO  TO  150 

1 AO  IF  ( NFL  AG  .EQ.  1)  GO  TO  AO 

IF  ( I OEL X  .LT.  0)  GO  TO  1A6 

RASTER  CO-OKO I N  ATE  OF  TIME  *  (ITIMEU)  -  TIMMIN)  /  (TIMMAX  -  T I  MM  IN)) 
*  (RASMAX  -  RASMIN)  ♦  RASMIN 

IA2  TIMOIF  =  TIMMAX  -  TIMMIN 
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00  1425  K 
ITIME(K)  = 
1425  CONTINUE 


C 


1«  NUMBER 
< (PLOTIM(K) 


TIMMIN1  /  TlMCIF)  *  1860  ♦  1020 


DO  144  J  =  1,  4 
C  DRAW,  TICK,  AND  LAPEL  AXES 

C 


143 


CALL  METAZ 
CALL  METAZ 
CALL  METAZ 
DO  143  K  * 
CALL  METAZ 
CCNTINUE 
CALL  METAZ 
CALL  METAZ 
CALL  METAZ 
CALL  MET AZ 


<9.  2,  I  AX ,  JAX •  KAX,  LAX) 
(9,7,KXT1CK,ITICK,KXTICK,JTICK) 

(9,9, IXTICK.KTICK, JXTICK.KTICK ) 

1,16 

(5,  0,  LOCX(K),  LOCY(K),  4,  LABEL ( K )  ) 

(5,  0,  IXPOS(l),  lYPOS(l),  3,  LABXAX) 

(50,  6,  1,  0,  0,  0) 

(5,  0,  I  XPOS(  2 )  ,  lYPOSm,  18,  LABYAXll.J)) 
(50,  6,  0,  0,  0,  0) 


RASTER  CO-ORDINATE  OF  CROSS  SECTION  =*  ((X(I)  -  PLMINX)  /  (PLMAXX  - 
PLMINX))  *  (RASMAX  -  RASMIN)  ♦  RASMIN 


DO  1435  K  =  1,  NUMBER 

JCROSS(K.J)  -  ( ( CROSEC ( K , J )  -  PLMINX)  /  RCSDIF)  *  2480  ♦  4C0 
1435  CCNTINUE 

DRAW  CURVE  AND  ADVANCE  FRAME 

CALL  METAZ  (8,  NUMBER,  ITIME(l),  JCROSS(l.J),  0,0) 

CALL  METAZ  (13,  1,  0,  0,  0,  0) 

144  CCNTINUE 

GO  TO  (  1443, 145) .NFLAG 

INITIALIZE  NEXT  PAIR  OF  FRAMES 

1443  XLEFT(l)  =  XLEFT ( NPLOT ) 

XR I TE (  1)  =  XRITE(NPLOT) 

YLEFT (  1)  =  YLEFT ( NPLOT  ) 

YR  I  TE (  1)  *  YRITE(NPLOT) 

PLOT  I M ( 1 )  =  PLOTIM< NPLOT) 

NPLOT  =  1 

1445  LABEL! !  )  =  TIMMAX 
DO  1447  K  *  2,  7 
LABEL ( K )  *  LABEL(K-l)  ♦  IDELT 
1447  CCNTINUE 

TIMMIN  =  LABEL!  1) 

TIMMAX  =  LABEL! 7) 

IF  ( PLOT  I M ( 1 )  .LE.  TIMMAX)  GO  TO  640 
630  IF  ( PLOT IM( 1 )  .LE.  TIMMAX  «■  DELSIX)  GO  TO  1445 
TIMMAX  =  TIMMAX  ♦  DELSIX 

GO  TO  630 

640  WRITE  (99,  6090)  ( LABEL ( K ) ,  K  *  1,  7) 

READ  (99,  6080)  ( LABEL ( K ) ,  K  -  1,  7) 

GO  TO  94 
145  CALI  PLTND 
C 

146  IF  ( NUMBUF  .GT.  1)  GO  TC  330 

147  IF ( I  OUT. EG. 8 )  CALL  EXIT 
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return 


150  IF  ( T AL2 ( NUM8UF  )  .LT.  TSTART)  GO  TC  AO 
IF  ( T AL2 ( NUMBUF  )  .LE.  TSTOP)  GO  TO  160 
MFLAG  =  2 
GO  TO  1  AO 
C 

160  NCHECK  =  A 

I T  AR2 ( NUMBUF )  =  -1 
NTGRN  =  ICAOATI  1  1 )  r  2 
IF  (NTGRN  .NE.  0)  GO  TO  162 
NTGRO  =  0 
NVHF2 ( NUMBUF  )  =  0 
NUHF2 ( NUMBUF  )  =  C 
NVHFO  =  0 
NUHFO  =  0 
GO  TO  290 
C 

162  NVHF  2 i NUMBUF  I  *  ICACATI21)  -  l 

NUHF2 ( NUMBUF )  =  NTGRN  -  NVHF 2  I NUMBUF ) 

IOENT  =  I  ACC  ♦  ICAOATI 17)  *  3 
NTCTE  =  ICADAT ( 8  ) 

DO  100  I  =  l.NTOTE 

I  DBASE  =  IOENT  +  (1-11*3 

K  =  IO0DAT (  1  )  +  1 

IF ( IDBCAT( 2 ) .EO.C  )  GO  TO  180 

NVTRG  =  I CAO A  T ( 2  1  ) 

IAMB(NTOTE)  =  I0BDAK3) 

JT  =  2 

I F ( IAMB(NTOTE).LT. NVTRG)  JT  *  1 
180  CONTINUE 
1B2  NTGRO  =  NTGRN 

NVHFO  =  NVHF  2 ( NUMBUF ) 

NUHFO  =  NUHF  2 ( NUMBUF ) 

1 8 A  IRORM  =  IRDRC2  ♦  IREL 

IF  (IRORM  .EC.  IREL)  GO  TO  290 
J  =  ICAOAT(IC) 

IF  (J  .EC.  0)  GO  TO  290 

GET  TIMES  1  RANGES  FROM  RORM 

I50MIL  =  I C AOA  T ( 1 3  I  *  3 
NRORM  =  ICAOATI  12)  *  N25MIL 
K25MIL  =  0 
00  200  I  =  1,  NRCRM 

CALL  RAOART  (FMRCRM,  IRORM,  IGMTH,  IGMTM,  GMTS,  TRTAL2II)) 
RANGE2 ( I )  =  ( RHDA  T ( 85  )  ♦  RD0ATC86) )  *  CCNRAN 
K25MIL  =  K25MIL  ♦  1 
IF  (K25MIL  .EC.  N25MIL)  GO  TO  190 
IRORM  =  IRORM  ♦  I25MIL 
190  K25MIL  =  0 

IR0RD2  =  IRCRO?  ♦  I50MIL 
IRORM  =  IRCRC2  ♦  IREL 
200  CONTINUE 
C 

C  DETERMINE  THE  BEST  MATCH  BETWEEN  MINOR  CYCLE  TIMES  AND  25  MS  DATA 
C  BLOCK  TIMES.  SAVE  THE  NUMBER  OF  THAT  MINOR  CYCLE  AND  THE  RANGE  FROM 
C  THAT  25  MS  DATA  BLOCK. 
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DTIME2  = 

icocoro. 

DEL  T l M  = 

(  ICADATI5I  /  1CCOCOO.  )  * 

ICADAT ( 7) 

DO  210  I 

*-  It  J 

TM I  NOR  = 

T AL2 ( NUMBUF )  ♦  ( I  -  1)  * 

CELT  IM 

DO  210  K 

=  lt  NRCRM 

DELTA  = 

DABStTMINOR  -  TRTAL2IK ) J 

IF  (DELTA  .GE.  DTIME2)  GO  TO  210 

DTIME2  = 

DELTA 

M l NCR  = 

I 

TRANGE  = 
CL  MTINUE 

R ANGE  ? ( K  ) 

GET  THE  RANGES  FCR  ALL  TARGETS  FROM  THE  M I  NCR  CYCLE  DETERMINED  ABOVE. 
DETERMINE  WHICH  RANGE  IS  CLOSEST  TO  THE  ONE  CHOSEN  ABOVE  AND  CALL 
THAT  THE  TRACKED  TARGET. 

ITRBAS  =  l ACC  ♦  ( I C ADA  T (  1 8 1  ♦  (MINOR  -  1)  *  (NTGRN  +2)  *  3)  *  3 
K  *  0 

DELRAN  =  ICOOOnO. 

IF  (IDENT  .NE.  I  ADD )  GO  TO  220 
GO  TO  (290,  21M.  KFLAG 
215  I  FLAG  =  2 
GO  TO  230 
220  I  FLAG  =  1 
KFLAG  =  2 
NTCT  =  NTCTE 

230  IF  (NTCT  .EC.  0)  GO  TO  290 
DO  270  1=1,  NTCT 
GO  TO  (250,  240),  I  FLAG 
240  K  =  K  +  1 
GO  TO  260 

250  I  DBASE  =  ICENT  ♦  (I  -  1 )  *  3 
J  =  I  DBCA  T ( 2 ) 

IF  (J  .EC.  0)  GO  TO  270 
K  =  K  +  1 

IAMB(K)  =  IDBDATI3)  *  3 
ITARG(K)  =  ICBDAT(l) 

IF  (  ITARG(K)  .OT.  1 )  GO  TO  255 
K  =  K  -  1 
GO  TO  270 

255  IPRI(K)  =  I  DO  D A  T ( 4  ) 

260  J  =  ITRBAS  ♦  lAMB(K) 

TDBRAN  =  IGET ( FMCTDB,  J,  1)  *  .0018737031 
DELTA  =  ABSITRANGE  -  TDBRAN) 

IF  (DELTA  .GE.  DELRAN)  GO  TO  270 
DELRAN  =  DELTA 
l  T  AR2 ( NUMBUF  )  =  ITARG(K) 

J  =  I  TAR  2 ( NUMBUF )  ♦  1 
I TYPE21 NUMBUF )  =  JTYPE(JT) 
l PR  102 ( NUMBUF )  =  IPRI(K) 

DR ANG2 ( NUMBUF  )  =  DELTA 
NCHECK  =  7 
270  CONTINUE 
NTOT  =  K 

290  DO  300  I  =  1,  NCHECK 

IF  (  I  BUFF ( NUMDUF , I )  .NE.  I BUFO (I)  )  GC  TO  310 
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300  CONTINUE 
GO  TO  40 

310  DO  320  1*1,  NCHECK 

IBUFOI  I )  =  I  BUFF ( NUMBUF , I  ) 

320  CONTINUE 

T  AG  T I  2 ( NUMBUF )  =  ( C  A  T  D AT  C  1 )  *  I2T22  ♦  CAT  DA  T ( 2  )  )  /  KCNTAG 

NUMBUF  =  NUMBUF  +  1 

IF  (NUMBUF  .LE.  MAXBUF)  GO  TO  40 

330  J  =  1 

340  IPAGE  =  IPAGE  ♦  1 

WRITE! I  OUT , 6  1 1 0  )  MESAGE, IPAGE 
WRITE! I  OUT , 60  30 ) 

00  370  I  =  1,  MAXLI2 

IF  (ITAR2IJ)  .NE.  -1)  GO  TO  350 

WRITE! I  OUT, 6040)  I  GMT H2 I J  ) , I GMTM<! t J ) , GMT S2 ( J ) , 

1 T  A  L  2  \  J  )  ,  TAG  T  I  2 ( J  ) ,  !IBUFF(J,K),  K  *  1,  4) 

GO  10  360 

350  WRITE! IOUT, 6040)  I GMTH2 ( J  ) , I GM TH2 ( J ) , GMT S2 ( J ) , 

1TAL2IJ),  TAGTI2U),  ( I  BUFF t J,K ) ,  K  *  1,  7),  DRANG2IJ) 

360  IF  (J  .GE.  NUMPUF  -  1 )  GO  TO  380 
J  *  J  ♦  1 

370  CONTINUE 
GO  TO  340 
C 

380  IF  (  ( NFL  AG  .EC.  2)  .AND.  (MFLAG  .EC.  2)  I  GO  TO  390 
NUMLIN  =  MAXLIN 
NUMBUF  =  1 
GO  TO  40 

390  IE! IOUT. EC. 8)  CALL  EXIT 
RETURN 
END 
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APPENDIX  E 

SUBROUTINE  CHEAD  PROGRAM  LISTING 


SUBROUTINE  CHEA0(*) 

C 

C  VERSION  04/01/71  R.H.  FRENCH 

C 

C  IHIS  SUBRUUlINE  IS  USfcD  !U  PRULfcSS  ALTA1R  CATALUb  TAPE  UA  f  A 

C  RECORDS. 

C 

C  IN  THIS  VERSION  FORMAT  TABLE  FMRR11  HAS  BEEN  ADDED.  THE  IDARR 

C  ARRAY  hAS  BEEN  PUT  IN  COMMON  SO  THAT  THE  RTS  AND  ARTP  VERSION 

C  OATES  C'.N  HE  POINTED  OUT  BY  THE  MAIN  PROGRAM.  THESE  OATE(S)  CAN 

C  ALSO  BE  USED  TO  TRIGGER  THE  FLOW  OF  THE  MAIN  PROGRAM. 

r 

INTEGER*2  ITEM 
INTEGER*2  IBU.TBU2 

COMMON/B'-AO/LN,  IFLG,  I AOO , FMRD 1 D , F MC A TF , FMC S AD , FMCMDB , FMCT I B , FMC I  OB 

1, FMCTDB,FMRORO,FMRORM,FMRDRT,FMGLOT,FMRRD5,FMAXSP,FMBIAS,FMR5CH 

2. FMRCHF, FMAACC.FMRR  1  1  .NAME ( 19),NK18),|X(18),I T6M( BODD ) 

COMMON/T I TLE/ IDARR (  10 ) 

DIMENSION  IBUF1 ( 2048) , I BUF 2 ( 204 B ) , MC A T F (18) .NAMEXI 19) 

EQUIVALENCE  ( F*RO I  0 , MC A TF ( 1 ) ) 

EQUIVALENCE  ( I PUF  1  (  1 )  ,  I BU) , C I BUF  2  I  1 ) , I BU2 ) 

DATA  NAMEX/'ROTO* , • C A TF • , ' C S AD ' , ’CMDB • , • C T I B • , • C  I DB  •  , » CTDB » , ' RDRD • 

1,  ' ROHM' ,' ROR T', 'GLGT'.'RR 05', ’AXSP', 'BIAS’,  ’R5CH’ , • RCHF ' , • AACC' 

2 ,  ' RR 1 1 ' , ' HERD ' / 

DATA  MAX/B192/, IFL/O/, IT/1/ 

NREC=0 
L  N=  1 

00  20  1=1,18 
MCATFI  I  )  =0 
NAME!  I  ) =NAMEX (  I  ) 

N I  I  I ) =0 

20  I  X  f I )  =  0 

NAME! 19)=NAMEX( 19) 

CALL  BREAOSILN,  IFjUFI,  1  BUF  2,  MAX,  IFL  ,  I  NDX  ,  L  EN  ,  I  FLG  ,  I  ACC) 

30  CALL  BREAD ( L  N  I 
NREC  =  NREC ♦  1 
I F (  IFLG.EQ.2  )GD  TO  55 
I F ( IFLG. EQ. 3. AND. I T.EQ.11G0  TO  30 
I F ( IFLG. EG. 3) GD  TO  90 
GO  TO  I  2  l , 22 ) , TNDX 

21  I F ( IBU/256.EQ. 1)G0  TO  90 
GO  TO  34 

22  I F ( IBU2/256.EQ. I )G0  TO  90 

34  CALL  NAMETI IAOD, NAMED) 

IF(NAMEO.EC.NAMEX(  19) )G0  TO  36 
WRITE(6,85INAMro 

85  FORMAT (  IX, A4,  •  FOUND') 

36  DO  IB  1  =  1,  19 

IF(NAME0.EQ.NAMEX( I ) )G0  TO  35 
18  CONTINUE 

GO  TO  30 

35  CALL  FORM! lAOO.ITEMI  IT/, IB. NAMED, NTEM,£70) 

I F ( NAMED  EC.NAMEX ( 19) ) GO  TO  £6 

DO  40  I  =  i,  18 

IF (NAME ( I ).EQ.NAMEO)CO  TO  60 
40  CONTINUE 
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GO  TO  30 

55  WR I TE ( 6  *  56 ) NREC 

56  FORMA T ( *  PARITY  ERROR  READING  FORMAT  RECORD*, 16) 

GO  TO  30 

60  MCATF (  I  )  =  IB 

WRITE(6,9)NAM£P,NTEM 

9  FORMAT  (  *  FORMAT  =  * , A4, '  STORED  IN  COMMON  NTEM*',I4J 
IX(  I  )  =  I T 
Nl ( I ) =NTEM 
I T  =  I T  *6*NTEM 

GO  TO  30 

66  CALL  HGRR! I  ADD, I TYP, 1DARR) 

WRITE(6.72)ITYP 

72  FORMAT,'  TYPE  *,I2) 

IF ( ITYP.EQ. I )GP  TO  30 
WRITE(6,68I 

68  FORMAT!'  TAPE  NOT  CATALOG  TAPE  JOB  TERMINATED  BY  CHEAD') 

RETURN  l 

70  WRITE(6*80)NAMCD,NTEM 

80  FORMAT ( '  NAME  =  ' A4 , '  NTEM  ='I5,»  *E  RRQR*  FORMAT  TABLE  LIMITED  T 
10  700  ITEMS  OR  FORMAT  TABLE  HAS  0  LENGTH') 

90  DO  81  1=1,  18 

IF ( MCATF ( I ).EQ.0)WRITE(6,82)NAMEX( I) 

82  FORMAT ( IX, AA,  •  NCT  FOUNO') 

81  CONTINUE 
WRITFI6,  150) 

150  FORMAT!'  lHEAO  COMPLETE') 

RETURN 

ENC 
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APPENDIX  G 

SUBROUTINE  CHEAD  OUTPUT 


I  D  =  CS 


type 

1 

LS  AD 

FUUNU 

FQHM4I 

=  CS  AO 

STORED 

IN 

COMMON 

NT  EM  = 

C  Tib 

FOUND 

E  UP  MAT 

=  CT  IB 

STORED 

IN 

COMMON 

NTEM  = 

UDB 

FOUND 

rUKMAT  =C IUb 

STORED 

IN 

COMMON 

NTEM  = 

LT  CB 

FOUND 

F  Ok  MA  T  =  C 1UH 

STOKtU 

IN 

COMMON 

NT  EM  = 

RURD 

FOUND 

FORMAT 

=  RDRD 

SI URED 

IN 

COMMON 

NTEM  = 

RURM 

FOUND 

FORMAT 

'  =  RDRM 

STORED 

IN 

COMMON 

NT£M  = 

RURT 

found 

FORMmI 

■=RDRT 

STORED 

IN 

COMMON 

NT  EM  = 

F  M  S  P 

FOUND 

PAC2 

FUUNU 

A  PC  3 

FOUND 

A  PCfa 

FOUND 

PAO  5 

FOUND 

PA  6  3 

FOUND 

1  RHD 

F  0  UN  D 

T  R  M  A 

FOUND 

1RTG 

FOUND 

TRMI 

FUUNU 

TRSP 

FOUND 

RDID 

FOUND 

FORMAT 

=  KD  I  D 

STORED 

IN 

COMMON 

NTEM  = 

AACC 

FOUND 

FORMAT 

=  A  ACC 

STORED 

IN 

COMMON 

NT  EM  = 

AM  P  1 

FOUND 

AMP2 

FOUND 

AMP  3 

FUUNU 

AMP4 

FOUND 

AM  P5 

FOUND 

AMP  6 

FOUND 

A  pc:  1 

FOUND 

AS  L  P 

FUUNU 

A  S  MP 

FUUNl) 

AT  RK 

r  UUNO 

AX  SP 

FOUND 

FORMAT 

=  AXS  P 

STOREO 

IN 

CUMMON 

NT  EM  = 

bCAL 

FOUND 

BIAS 

FOUND 

FORMAT 

=  61  AS 

ST  ORED 

IN 

COMMON 

NT  EM  = 

bSMC 

FOUND 

CATF 

FOUND 

FORMAT 

=  C  AT  F 

STORED 

IN 

CUMMUN 

NT  EM  = 

CH  AF 

found 

CM  Ob 

found 

F  URMAT 

=  CMDll 

ST URED 

IN 

COMMON 

NT  EM  = 

GLCT 

FOUND 

FORMAT 

=  GLO  T 

STORED 

IN 

COMMON 

NTEM  = 

NO  Ml 

FUUNU 

NO  M2 

FOUND 

N  JM 3 

FOUND 

N0M4 

FUUNU 

btiJT 

FOUND 

PHA1 

FOUND 

PHA2 

FOUND 

PHA3 

FUUND 

PHA4 

FOUND 

6 

PH  A  5 

FOUND 

PHA6 

FOUND 

3 

KCHF 

FOUND 

FORMAT 

'RCHF 

STORED 

IN 

COMMON 

5 

R4CH 

FOUND 

R5CH 

FOUND 

4 

FORMAT 

=R5CH 

STORED 

IN 

COMMON 

RR04 

FOUND 

18 

RRCI5 

FOUND 

FORMAT 

*RR  05 

STORED 

IN 

COMMON 

i2  7 

RR06 

FOUND 

RRU7 

FOUND 

10 

RRCfl 

FUUND 

RR09 

FOUND 

ERRO 

FOUND 

SCAN 

FOUND 

X  S  EC 

FOUND 

DRG1 

FOUND 

DR  02 

FOUND 

URG3 

FUUND 

SA  SC 

FOUND 

KR1 1  NDT  FOUND 

CHEAD 

COMPLETE 

4 

8 


29 

19 

23 

51 

7 


NT  EM  =  6 

NT  EM  =  32 

NTEM  =  8 
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APPENDIX  H 

SUBROUTINE  LTIME  PROGRAM  LISTING 


SUBROUTINE  LTIMEILCT,IH,IM,IS,IT) 
DIMENSION  LOT H ) - ICON! 5 ) 

DATA  I  CON/ Z  1C, 7 1C0,Z 1 000 , Z 10000, Z 1 COOOC/ 
IS1  =  MODILOTI?), ICON!  1)) 

I  S  10  =  MCDILCT  I  2  )  /  ICON  (  1),  ICON!!  1 )  ) 

IS  =  IS10*10  +  ISi 

IM1  =  MODILOTI ?)/ICCN(2),  ICON! 1) ) 

iMIO  =  MOCILCTl 2) / I  CON ( 3)  ,  ICON! 1 ) ) 

IM  =  I M  1C*  1C  <■  I M  1 
I H  1  =  MODILOTI  ?)/ICCNU)  ,  ICONI  D) 

I H 10  =  MCDILCT I  2 ) / ICON  I  5 ) , ICONI  1  )  i 

IH  *  I H 10  *  10  +  I H 1 

ITU  =  MODILOTI 1) , ICONI  1  )  ) 

ITT  *  MODILOTITJ/ICONI  1), ICONI 1) ) 

ITH  =  MOOILOT I  1 ) / ICON ( 2 ) « ICON( 1 ) ) 

IT  -  ITU  +  ITT* 10  ♦  ITH* 100 

return 

END 
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APPENDIX  J 

SUBROUTINE  GMTTAL  PROGRAM  LISTING 


SUBROUTINE  GMTTAL  (IGMTH,  IGMTM,  GMT  S  t  TALI 

COMMON  /LAUNCH/  TLCNCH 

INTEGER  *  2  IGMTH,  IGMTM 

REAL  *  8  GMTS,  TAL 

REAL  *  8  TLONCH 

DATA  I2T20  /10AP576/,  I2T16  Z65536/,  I2TS2  /A096/,  I2T8  /256/, 
1 1  2  T  A  /  1 6  / 

COMPUTE  SECONOS  AFTER  LAUNCH 

100  TAL  =  IGMTH  *  3600  ♦  IGMTM  *  60  ♦  GMTS  -  TLCNCH 

CHECK  TO  SEE  IF  TAL  IS  NEGATIVE.  IF  SO,  CHECK  FCR  MIONITE  LAUNCH 

IF  (TAL. LT. 0.1  GC  TC  10 
RETURN 

10  IF< CABS! TAL 1 .GF.837C0.0)  GO  TO  11 

TAL  =  IGMTH  ♦  36C0  ♦  IGMTM  *  60  ♦  GMTS  -  TLCNCH 
RETURN 

11  TAL=TAL*86A00. 

RETURN 


ENTRY  CATIME  (FMCATF,  IAOD,  IGMTH,  IGMTM,  GMTS,  TAL) 

J  =  I G  E  T ( FMCATF ,  IAOO,  3) 

IGMTH  =  ( J  /  I2T20)  ♦  10  ♦  (J  -  (J  /  I2T20)  *  I2T20)  /  I2T16 

IGMTM  *  ((J  -  (J  /  1 2  T 1 6 1  ♦  I2T16)  /  I2T12)  *  10  ♦ 

1 ( J  -  (J  /  I  2  T  1  2 1  ♦  I  2  T  1 2  I  /  I2T8 
GMTS  =  ((J  -(J  /  I2T8I  *  1 2T8 1  /  I2TA)  *  10  ♦ 

1 ( J  -  (J  /  1 2T A 1  *  I2TAI 
J  =  I  GET ( FMCATF,  IACO,  A) 

GMTS  =  GMTS  ♦  (J  /  I2T8I  *  .1  ♦  ((J  -  (J  /  I2T8)  *  I2T8)  /  I2TA) 
1*  .01  ♦  (J  -  (J  /  1 2T A  J  *  I2TA)  *  .001 

CONVERT 

GO  TO  IOC 


ENTRY  RADART  (FMRORM,  IKCRM,  IGMTH,  IGMTM,  GMTS,  TAL) 

J  =  I GE  T ( FMRORM ,  IRCRM,  l) 

IGMTH  =  (J  /  17T12I  *  10  ♦  (J  -  (J  /  I2T12)  *  I2T12)  /  I2T8 

IGMTM  =  ((J  -  (J  /  I 2T8 )  *  I2T8)  /  I2TA)  *  10  ♦ 

1(J  -  ( J  /  I 2TA )  ♦  I  2 T  A  I 
J  =  I GE  T ( FMRORM ,  IRDRM,  21 

GMTS  =  ( J  /  1 2  T 1 6 )  *  1C  ♦  (J  -  (J  /  I2T16I  *  I2T16)  /  I2T12  ♦ 
1((J  -  (J  /  I2T12)  *  I2T12)  /  I2T8)  *  .1  ♦ 

2 ( ( J  -  (J  /  I2TP)  *  I2T8)  /  1 2  T  A )  *  .01  ♦ 

3 ( J  -  (J  /  1 2T A )  *  I2TA)  *  .001 

CONVERT 

GO  TO  100 
END 
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APPENDIX  K 

SUBROUTINE  REFC  PROGRAM  LISTING 


* 


t 


100 

120 

200 

220 


300 

320 

340 

400 


SUBROUTINE  R E F C ( E , R » 0 E E , ORR  1  ypncinM.  6/16/7H 

DIMENSION  CE(16,8),0R(  16,8)»E0(L6)tRD(8)  VERSI0N-  6/16/70 
OA  T  A  DE / 0. 0  ,0. 0  ,0.0  ,0.3  ,0.0  ,0.0  ,0.0  ,0.0  i 

10.0  ,0.0  ,0.0  ,0.0  ,0.0  ,0.0  ,0.0  ,0.0  ,0.0313,' 
20. 0303,0. 0292, 0.02e7, 0.0282, 0.0272, 0.0 262, 0.0253, 0.02 43,0.0223, 
30.0214,0.0195,0.0171,0.0135,0.0075,0.0  ,0.0937,0.0848,0.0770, 

40.0732,0.0694,0.0627,0.0571,0.0522,0.0480,0.0412,0.0385,0.0337, 
50.0278,0.0205,0.0105,0.0  ,0.1850,0.1520,0.1250,0.1140,0. 1050 , 

60.C904, 0.0 795, 0.0 708, 0.0636, 0.05 23, 0.0478, 0.0405, 0.0323,0.0 229, 
70.0114,0.0  ,0.5310,0.3070,0.2120,0. 1830,0.1600,0.1280,0. 1060 , 

80.0899,0.0780,0.0612, 0.0550,0.0455,0.0354,0.0246,0.0120,0.0  , 

90.7550,0.3720,0.2400,0.2020,0.1750,0.1370,0.1120,0.0942,0.0811 , 
AO. 06 31, 0.0566, 0.0466, 0.036 1,0. 0250, 0.0122, 0.0  ,0.9120,0.4110, 

BO. 2560, 0. 2 14C , n . 1840,0. 1420, 0.1 150, 0.0967, 0.0830, 0.064 3, 0.0575, 
CO. 0472, 0.0 365,0. 0252, 0.0122, 0.0  ,0.9700,0.4200,0.2600,0.2200, 

DO. 1900,0. 1460,0. 1170,0.0980,0.0840,0.0653,0.0584,0.0478,0.0369, 
E0.0254, 0.0123,0.0  / 

OATA  DR/  0«0*  0«0f  0«0*  0*0*  0«0»  0*0«  0«0t  0«0«  0«0« 


0.0 


xs.uv  0.0,  0.0,  0.0,  0.0,  0.0,  C.O,  22.6,  21.5,  20.4,  19.9, 

19.4,  18.5,  17.6,  16.8,  16.1,  14.8,  14.2,  13.2,  12.0,  10.4,  8.6, 

7.7,  67.3,  57.9,  50.2,  47.0,  44.1,  39.3,  35.4,  32.1,  29.3,  24.8, 

22.9,  19.7,  16.3,  12.7,  9.4,  8.1,132.0,  98.5,  77.4,  69.7,  63.2, 

C  -I  J  i  t  ■>  n  i  t  t  ~  «  ... 


8  •  6  I 


1 

2 
3 

n  «-«-•*»  *  '•  '  »  1LI  II  »  Oi  If  IJ^tUf  70*  3v  (  r«S| 

5  52.9,  44.7,  38.4,  33.4,  26.4,  23.9,  20.1,  16.4,  12.7,  9.4,  8.1, 

6340.0,  167.0, 10 ’.0,  86.  1,  73.4,  56.7,  46.2,  38.9,  33.6,  26.4,  24.0, 

7  20.2,  16.4,  12.8,  9.5,  8.2,405.0,170.0,104.0,  86.3,  73.6,  56.8, 

8  46.3,  38.9,  37.7,  26.5,  24.1,  20.3,  16.5,  12.8,  9.5,  8.2,421.0, 

9171.0,104.0,  86.6,  73.9,  57.1,  46.4,  39.0,  33.8,  26.8,  24.3,  20.5, 

A  16.6,  13.0,  9.8,  8.4,446.0,172.0,105.0,  87.4,  74.0,  58.0,  46.6, 

8  39.2,  34.0,  27.0,  24.6,  20.7,  16.7,  13.0,  10.0,  8.4/ 

DATA  E0.RT0EG/0. 01, 2. 0,4, 0,5. 0,6. 0,8. 0,10. 0,12.0, 14.0, 18., 20., 

124. ,  30.  ,40. ,60. ,90. ,5  7.295  7  8/ 

DATA  R 0/0.01 , 10., 30. ,60. ,200. ,400. , 1 000. , 2000 . / 

IFIR.LE.O.CIGO  TC  300' 

RG  =R /60  80 . 2  7 
OD  100  IEO-2,15 
I  =  1 7- I FO 

IF ( E.GE.EOI I ) ) GO  TC  120 
CCNT INUE 
1=1 


00  200  JR0  =  2 , 8 
J= 10-JRD 

IFIRG.CE.ROI J))GC  TC  220 

CONTINUE 

J=1 

IFIJ.EC.BIGC  TC  34C 

ZR=ALOGI RG/ROI J) ) / 4LOG ( RD 1 J+ 1 ) /RD ( J ) ) 

ifie.le.o.oigd  re  320 

ZE=ALCG( E/EDI  I  I ) /ALCG ( ED ( 1+ l ) /ED (1)  ) 

OEl*MCE(  !♦!  J)-DE(  I ,  J  )  )*(  l.-ZR)4(0E(  1  , J+ 1 ) -DE (  I ,  J  )  ) *Z  R ) *Z  E 
0E2=( (0E(  I, J»! I-OEU, J) )*» l.-ZE )*(DE U  +  l, J+  1)-DE( I, J+l ) )*ZE)*ZR 
0EE=0El+0E2»CEI I, J) 

OR l = I (DR( 1+ 1 , J ) -CR ( I,J) 1  *  ( l.-ZR)+(DR( I ,J+1)-DR(I,J))*ZR)*ZE 
OR  2=  (  (  CR  (  I ,  J*  D  -CR  (  I ,  J  1  )  *  (  l  .-ZE  )♦  ( DR  (  I  + 1 ,  J+  l  l-DR  I  I ,  J+  1  ) )  *2  E)  »ZR 
ORR  =  ( DR  1  +  CR2+0R I  I , J  )) 

GO  TC  400 
DE  E  =  0 . 0 


0RR=0. 0 
GO  TO  400 

OEE=OEI1, JI+IOEI I,J+l)-OE< 1,J) 1*ZR 
ORR=OR ( I , J ) ♦ (OR ( I, J+l) -OR! I , J  » )*ZR 
GO  TC  400 

OEL  T= ( E-ED ( I ) l/IEDI I  +  l ) -E D (  1  )  ) 
OEE=l>a.T*(ClF(  I+1,J)-DE(  I,  J)  )+DEU  ,  J) 
DRR=0EIT8(DRI !+l, J )-DR( I, J))+ORU , J) 
RETURN 
ENO 
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APPENDIX  L 

SUBROUTINE  BZERO  PROGRAM  LISTING 


JU,  J 


START 

0 

ENTRY 

BZERO 

OSINO 

*,  15 

B7ERC 

B 

H  5 

DC 

X*  05*  ,CL5  *  BZEFC* 

H6 

ST1 

14,  12,  12  (1?) 

I. 

6,  0(1) 

NVC 

WORD  ,0  (6) 

L  E 

4,  WORD 

tPER 

6,4 

CE 

6  ,  =E  1.  ' 

BL 

R2  0 

AU 

4,  X  6 

STE 

4 , WORD 

I, 

7  ,  WOR  D 

SR 

9,4 

LTR 

7,7 

SLL 

7,9 

BNL 

HI  0 

LA 

9,  1 

HI  0 

SRL 

7,9 

SR 

2,2 

AH 

2,  -X '  010  3' 

SRDI, 

2,2 

SRL 

3,30 

SLL 

7,0(3) 

SFL 

7,2 

SRDL 

2  ,8 

OR 

7,3 

SLL 

9,31 

OR 

n  ,9 

SER 

2,2 

ST 

7, WORD 

AE 

2 , WORD 

STE 

2,  WORD 

WVC 

0  (4  ,6  )  ,  WORD 

H20 

RETURN 

(14,  12) 

WORD 

DS 

IE 

X6 

DC 

X*  46  000000  • 

END 
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APPENDIX  M 

FUNCTION  IBIT  PROGRAM  LISTING 


FUNCTION  IBIT  (ITEM,  NBIT) 
INTEGER  *  2  I TFM 
LIMIT  =  NBIT  -  2 
J  =  0 

IOO  IF  (ITEM  .EQ.  (2**J)  )  GC  TO  120 
IF  (J  .GE.  LIMIT)  GO  TO  110 
J  =  J  ♦  1 
GO  TO  100 

110  IBIT  =  NBIT 
RE  TORN 

120  IBIT  =  J  +  1 
RETURN 
END 


